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This paper is concerned with an analytical model of traffic flow on networks consistent 
with the Lighthill-Whitham- Richards model [39J [47] . The model captures the formation, 
propagation and dissipation of physical queues as well as spillbacks on a network level. 
Similar models have been mainly studied either in discrete-time p~2j [13] [51], or at the 
level of specific junctions [29] [30j [34]. In this paper, we take one major step towards an 
analytical framework for analyzing the network model from two approaches: (1) we de- 
rive a system of differential algebraic equations (DAEs) for the explicit representation of 
network dynamics. The DAE system is the continuous-time counterpart of the link trans- 
mission model proposed in [51] ; (2) we conduct mathematical analysis of the continuous- 
time model and discuss solution existence, uniqueness as well as continuous dependence 
on initial/boundary conditions. The proposed work has a positive impact on analytical 
Dynamic Traffic Assignment models in that it not only provides theoretical guarantee 
for continuous-time models, it also brings insights to the analytical properties of discrete 
models such as the cell transmission model p~2j [32] and the link transmission model [51] . 

1 Introduction 

This paper focuses on the Lighthill-Witham- Richards model [39], [47] on a vehicular network. 
We model the traffic dynamics on a link level with the following first order hyperbolic partial 
differential equation (PDE) which describes the spatial-temporal evolution of density and flow: 



where p(i, x) : [0, +oc) x [a, b] — > [0, pj] denotes density, f(p) : [0, fP am ] — > [0, C] denotes 
flow. fP am is the jam density, corresponding to a bumper-to-bumper situation, C is the flow 
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capacity. The function /(•) describes explicitly the density- flow relation and is referred to as 
the fundamental diagram. 



Classical mathematical results on the first-order hyperbolic equations of the form (1.1) 
can be found in [3J. For a detailed discussion of numerical schemes of conservation laws, we 
refer the readers to [241 EBJ. A well-known discrete version of the LWR model - the Cell 
Transmission Model (CTM) - was described in jT2l [13]. The PDE-based models have been 
studied extensively also in the context of vehicular networks, with a selected list of references 

including ©Eiiainiiiainai^ 

1.1 Network models 

The fluid-based models on networks turn out to be much more complicated and subtle than a 
single conservation law. The drivers' behavior around merging and diverging intersections will 
affect the boundary data on each link incident to that intersection in a highly nonlinear way 
[22J . Let us consider a network consisting of just one junction with 771 incoming links I\ , . . . , / m 
and n outgoing links I m +i, . . . , Im+n- The notion of a weak solution at this junction associated 
with the system of conservation laws [27] is analogous to the one for a single conservation law. 
That is, the densities on each link k — 1, . . . , m + n satisfy 

*2±? r+oo r b k d fa d(j)k 

Pk^rr + fk{Pk)^—dxdt = (1.2) 
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for all smooth 0& with compact support in (0, +oc) x (a&, k — 1, . . . , m + n, that are also 
smooth across the junction: 

</>i(tbi) = <j>j(t,a,j), -7j^(*> b i) = "^"(*' for a11 * 

where i = l,...,m, j = m + 1, . . . , m + n. For any weak solution Pk(t, •)? ^ = 1? • • • j m + n 
with bounded variations for each t, the Rankine-Hugoniot condition holds: 

m m+n 
i=l j=m-\-l 



Notice that (1.3) describes the conservation of flow through junction. However, the solution at 



a junction remains largely underdetermined with (1.2) and (1.3). In order to isolate a unique 
solution, various entropy conditions were proposed and analyzed, see pTU| fT3| l27| l30| H3] . 

The analytical junction models are best analyzed using the Riemann problems. A Riemann 
problem at a junction is an initial value problem corresponding to an initial data that is 
constant on each link. The Riemann problem and its solution method, the Riemann Solver, are 
extremely helpful to construct and analyze junction models in continuous time. The Riemann 
Solver associates, to each set of Riemann initial conditions, a set of boundary conditions such 
that the solution at the junction is obtained by solving a set of initial/boundary value problems 
on each link connected by the junction, (see jTBJ E2| for an overview of the Riemann Solvers). 
Various Riemann Solvers for merge and diverge junctions associated with different rules or 
policies, have been derived in the literature [TOj, [131 HH1 EH [29l [30] . The Riemann Solver can be 
explained easily with the notion of demand-supply put forward in [3H [35l [36]. The demand 
(supply) for an incoming (outgoing) link is a function of the density p(t, fy— ) (p(t, %•+)), 
which measures the maximum flux that an incoming (outgoing) link can send (receive). 
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1.2 Analytical models 

Despite the fast growing literature on junction models of various types, the analytical proper- 
ties of the system of conservation laws coupled with junction models are less understood. Solu- 
tion properties such as existence, uniqueness and continuous dependence on initial/boundary 
data have not yet been confirmed by discrete-time models. [IB] found out that capturing 
realistic traffic dynamics, such as queue spillback, may violate the requirement on solution 
existence. In [35] , it was shown that under congested conditions, route travel times obtained 
with the CTM may become discontinuous. In addition, analytical representation of traffic flow 
that adequately replicate traffic theoretic relationships and yield well-behaved mathematical 
formulations are currently unavailable [2]. 

The scope of this paper is the continuous-time fluid-based model on vehicular networks. 
We aim towards an analytical representation of the dynamics on the links and through the 
junctions, while providing theoretical guarantee of solution existence, uniqueness and continu- 
ous dependence on initial/boundary value conditions. Such topics are crucial to a well-defined 
and well-posed theory of fluid-based traffic flow models. They also serve as the very founda- 
tion of analytical dynamic traffic assignment (DTA) models such as dynamic user equilibrium 
(DUE) and dynamic system optimum (DSO). The DUE is formulated as either a variational 
inequality [19], or a nonlinear complementarity problem [25, 50J or as a differential variational 
inequality [201ET]. The DSO problem is usually formulated in continuous-time, as an optimal 
control problem [18] . or a variational inequality [46] . Solution existence of all the aforemen- 
tioned formulations typically requires the continuity of the delay operator arising from the 
dynamic network loading sub-problem [19J. The continuous dependence on initial/boundary 
data of network solutions as we put forward in this paper plays a key role in showing continuity 
of the effective delay operator in a Hilbert space. 

Our main technique for showing solution existence and continuous dependence is the use 
of a wave-front tracking algorithm and suitable approximations and functionals to control 
the total variation of the fluxes across the junction. The wave- front tracking algorithm is a 
basic technique for constructing solutions to system of conservation laws, based on piecewise 
afflne flux function and piecewise constant initial/boundary datum. We refer the readers to [3J 
and [28J for a general theory of conservation laws and for a discussion of wave-front tracking 
algorithms. The main difficulty in constructing solutions to system of conservation laws via 
wave-front tracking arises from the control of total variation of the approximate solutions 
Pv(-, •)> v > 1. For a single conservation law, the total variation of the solution p v {t, •) is 
decreasing with respect to time, and the number of wave interactions is finite, see [3J. While 
for systems of conservation laws which is the case of network flow models, the construction 
of wave-front tracking approximation is more subtle due to the interaction of waves on an 
link with junctions. In particular the bounds on number of waves, interactions and total 
variations are much more difficult to estimate. In [10J, an example was shown where a single 
wave of arbitrarily small variation generates, after interacting with the junction, waves whose 
variations are bounded away from zero. Such arbitrarily large magnification of total variation 
is unique to conservation law systems. As a consequence, one could no longer expect any 
bound on the total variation of the solution in terms of the total variation of the initial data. 
In order to proceed, we need to define carefully the Riemann Solver we want to employ in the 
construction of a wave-front tracking solution. 

Another contribution made in this paper is an analytical representation of the dynamics on 
a network- wide level. In general, the LWR-based models do not admit closed form solutions 
even for a single conservation law. Nevertheless, by assuming a triangular fundamental dia- 
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gram and an initially empty network, we are able to implicitly solve the dynamics on each link 
given the boundary data using a variational argument. Such type of argument has been used 
for the analysis of traffic flows on a single (possibly inhomogeneous) link by Daganzo |14( [T5] 
and Newell [42J. It has also been applied to traffic research extensively [4j |5j [8j |9l [43J [51] . In 
our formulation, we rely on the analysis of a separating shock wave on each link [6j. Given 
that the link is initially empty, the separating shock will be the unique shock wave on the link 
that separates the free flow phase and the congested phase. Depending on the left and right 
boundary conditions, the separating shock can reach and stay in the entrance or the exit of 
the link for some time period. When this happens, we call the separating shock a latent shock. 
The case where the separating shock reaches the entrance of a link is identified as spillback 
since the current link is dominated by congestion and as a result, the supply at the entrance is 
reduced, such reduction of supply poses further constraints on the junction solution and affects 
the upstream links. We propose an analytical way of detecting the latent shock based only on 
the cumulative arriving and departure curves at the two boundaries of the link. Finally, given 
appropriate Riemann Solver, the dynamics on a network level can be described explicitly as a 
system of differential algebraic equations (DAEs). The solution properties established in this 
paper will apply to such DAE system, as it is just a special case of what we will analyze using 
wave-front tracking algorithm in this paper. 

We make an important note that the resulting DAE system is the continuous-time coun- 
terpart of the link transmission model (LTM), which was originally proposed as a numerical 
scheme in [51]. In this study, we do not attempt to develop new traffic flow models or com- 
putational methods. Rather, we intend to establish a solid theoretical background for the 
continuous-time traffic flow model, which is the foundation of analytical DTA models. The 
properties of continuous models established here will also provide insights on discrete-time 
models such as CTM and LTM. The analytical formulation of network dynamics as well as 
rigorous theoretical results reported in this paper constitute a major step towards a mathe- 
matically sound theory of network flows. 

1.3 Organization 

The rest of the paper is organized as follows. In Section [2j we will present an analytical 
framework for analyzing the continuous-time model on networks. In particular, we rely on a 
continuous variable indicating flow and a binary variable indicating the regime. The system 
can be greatly simplified by analyzing the separating shock [6J on each link, and their interac- 
tions through the junctions. Section [3] summarizes the discussion in Section [2] and present an 
explicit DAE system for describing and predicting the network dynamics, including spillback. 
The DAE system is the continuous-time counterpart of the Link Transmission Model [51J. Sec- 
tion [4] is devoted to establishing analytical properties of the resulting continuous-time system. 
The argument we use is generalizable to a system with any BV (bounded variation) initial 
and boundary conditions. In particular, we show solution existence, uniqueness and continu- 
ous dependence on initial/boundary data, in the case of two Riemann Solvers introduced in 
Section |2| In Section [5| we will present a numerical example of the model on a network. Note 
that we do not attempt to replicate the LTM algorithm in this numerical study. Rather, the 
point is to show the features of the proposed continuous-time model such as separating shock, 
latent shock as well as spillback. 
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2 Formulation of the Link-based Kinematic Wave Model on 
networks 

In this section, we will present a continuous-time fluid-based model for a road network, as- 
suming the triangular fundamental diagram. This model relies on a new set of state variables, 
and the dynamics on a link level can be inferred from the boundary conditions via Lax-Hopf 
formula [HI [TU [T71 [32] . Note that the time-discretization of such system is precisely the 
link transmission model (LTM) proposed in [5T] , Unlike the cell-based models [El EH], the 
Link Transmission Model does not require explicit computations in the interior of the links, 
rather, it relies on the variational argument proposed by Newell [42J to determine the send- 
ing/receiving flows on a link level. The continuous-time formulation proposed in this paper 
complements the LTM and provides further insights of the LTM in the context of analytical 
models. In particular, we will establish properties of continuous-time solution to this system 
later in Section HJ 

2.1 State variables of the system 

Consider a link represented by an interval [a, 6], with b — a = L > 0. In the formulation of 
LKWM model, we will choose flow q(t, x) and regime r(t, x) as the state variables for the 
system. It is obvious from Figure [T] that a single value of q corresponds to two traffic states: 
1) free flow phase (r = 0) and 2) congested phase (r = 1). Hence the pair (q(t, x), r(t, x)) E 
[0, C] x {0, 1} determines a unique traffic state. The following map is then well-defined 

</>(•) : [0, C] x {0, 1} [0, p> am ], (q, r) H- p (2.4) 
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Figure 1: Illustration of the new state variable (</, r). 
2.2 Traffic flow through junctions 

Various junctions models associated with the LWR model have been developed in different 
contexts. For instance, in one research line, [27J approached the junction model by invoking the 
concept of weak solution for a system of conservation laws across the junction, as well as the 
notion of entropy condition which amounts to maximizing a certain quantity over all admissible 
boundary fluxes. [10] solved a Riemann problem on a junction using the concept of Riemann 
Solver, i.e. a map associating, to each set of Riemann initial conditions, a set of boundary 
conditions such that the solution at the junction is obtained by solving initial/boundary value 
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problems on each link connected to the junction. This approach was extended in subsequent 
papers, see [Tf \16\ HO]. 

In another research line, the junction model can be described less technically using the 
demand-supply functions [35], [36] . The demand function determines, based on the traffic con- 
ditions immediately upstream to the junction, how much traffic can leave the upstream links; 
while the supply function determines, based on the traffic conditions immediately downstream 
to the junction, how much traffic can enter the downstream links. Various rules and policies 
can be associated with the demand and supply functions to account for different modeling as- 
sumptions. A selected list of different merging and diverging models based on demand-supply 
functions include [HI Eg EQl EH S3] 

In this section, we will present the Riemann Solvers [10J for merging and diverging junctions 
of simple types. In order to obtain an analytical representation of the junction models, we 
consider only two types of junctions depicted in Figure |2| More complicated junctions can be 
decomposed into a set of these two basic ones by introducing virtual links. Similar techniques 
have been applied in [T3] . The Riemann Solver proposed in [IQ] must satisfy the following 



assumptions. 

(Al) There exists a given distribution matrix A that prescribes the turning percentage of 
drivers from incoming links to the outgoing links, when the number of the outgoing 
links is > 2. 

(A2) Subject to (Al), the flux is maximized through the junction. 

Remark 2.1. As discussed in [10], the above conditions satisfied by a Riemann Solve amount 
to a linear programming problem, where the flux is maximized with constraints given by the 
demand and supply of incoming and outgoing links, as well as the distribution matrix A. 

2.2.1 Diverge model 

Consider the diverging node of Figure [2| with one incoming links I\ and two outgoing links 
I2, ^3- Each link 1{ is characterized by a spatial interval [a^, 6j], a fundament diagram fi(-) 
and a flow capacity C{. The demand function D\{t) and the supply functions Si(t) for i = 2, 3 
are defined as follows 



Ii 

> 




Figure 2: Basic diverging (left) and merging (right) nodes 




if n(t, 61-) = 1 
if n(t, 61-) = 



(2.5) 




if n(t, a,i+) = 1 

if ri(t, a,+) = 



(2.6) 
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Recall that the regime variable satisfies r = 1 if the traffic is in congested phase and r = 
if traffic is in free flow phase. 

The diverging model requires the knowledge of a (possibly time-dependent) vehicle turning 
rate E i — 2, 3 such that + Q?i,3 = 1. These parameters can be given either 
endogenously or exogenously, depending on the context in which the model is considered. 

We employ the diverging model satisfying first-in- first-out, proposed in [101 U3] , the flow 
leaving the incoming link q ou t,i(t) can be easily solved by the LP and is given by 

= min\ Dl (t), (2.7) 

The inflow on link U is given by qi n ,i(t) = qlij • q ou t,l(t), i = 2, 3. 



2.2.2 Merge model 

Consider the merging node of Figure [2j with two incoming links 1^ I5 and one outgoing link 
Iq. The demand and supply functions for incoming and outgoing links respectively, are defined 



in the same way as in (2.5), (2.6). It has been pointed out in [10] that for this merging junction 
(more generally, for junctions with more incoming links than outgoing links), the LP problem 
defined above does not have a unique solution. Therefore, we need to introduce a right-of-way 
parameter p E (0, 1) and a priority rule in the case of two incoming links ^4, I5 and one 
outgoing link Iq. 

(Rl) The actual fluxes q ou t^ q ou t,5 out of links ^4, I5 satisfy 

P ' qout.A = qout,5 

for some right-of-way parameter p E (0, 1). 

However, one should notice that the rule (Rl) is sometimes incompatible with the assumption 
(A2). See Figure [3] for an illustration. 



^out,5 





^out^ 



Figure 3: Left: Rule (Rl) is compatible with (A2), therefore there exists a unique point 
satisfying both. Right: Rule (Rl) is incompatible with (A2), in this case, the Riemann Solver 
chooses the point that is closest to the line with slope p among all points that maximizes the 
flux. 

Whenever there is a conflict between (Rl) and (A2), we will respect (A2) and relax (Rl) 
so that the solution to the LP problem is chosen to be the point that approximates best the 
rule Rl among the set of points that attain the maximal flux. The Riemann Solver defined 
in this way uniquely determines a solution to the Riemann problem at junctions of merging 
type. 
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2.3 Shock formation and propagation within the link 

In the previous section, we have expressed the merging and diverging models in terms of 
demand and supply, which are functions of the regime variables r at the upstream and down- 
stream boundaries. Thus the last missing component of a complete network model is the 
determination of the traffic states (free flow or congested) on the boundaries of each link. 
This is the primary goal of this section. 

The key point of our analysis will be the location of a separating shock within the link. 
The separating shock is a generalized characteristic (for a precise definition and theoretical 
results, we refer the readers to which separates each link into two zones of free flow 

traffic (r = 0) and congested traffic (r = 1). 

We begin with showing that if the network is initially empty, then there can be at most 
one separating shock within each link 

Lemma 2.2. For every link I{ and any solution yi(t, x) — x), r^(t, xj) with yi(0, x) — 

(0, 0), the following statement holds: 

1. For every t > 0, there exists at most one x*(t) E (a^, b{) such that n(t, x*(t)—) < 
ri(t,x*(t)+) 

2. For all x E [a^, bj\, 

ri(t, x) = 0, if x < x*(i) 
Ti(t, x) = 1, if X > x*(t) 



Proof. See [6J. 



□ 



Remark 2.3. Lemma holds also for the case of network completely in free regime initially. 
In other words it is sufficient to ask for r^(0, x) = for every i and x. 



By Lemma pO} the separating shock will emerge from the downstream boundary bi of the 
link, and propagate towards the interior of the link. The speed of this separating shock is 
given by the Rankine-Hugoniot condition [T7] , see Figure |4j 
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Figure 4: Example of the separating shock. Left figure: the t — x domain is separated by the 
shock into two phases. Right figure: the speed of the separating shock ^x*(t) is given by the 
Rankine-Hugoniot condition, i.e., the slope of line segment on the fundamental digram that 
connects (qi e ft, 0) and (q rig hu !)■ 
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It is clear that as long as the separating shock remains in the interior of the link 1^ the two 
boundary conditions do not interact. As a result, the exit of an incoming link remains in the 
congested phase, while the entrance of an outgoing link remains in the free flow phase. Con- 
sequently, the demand-supply functions and hence the junction models are expressed entirely 
with endogenous parameters C^, ctij. On the other hand, if the separating shock reaches the 
boundary, it becomes a latent shock. We distinguish between two extreme cases. 

i) The the shock reaches the exit, i.e. x*(t) = b{. In this case, the current link is entirely 
in free flow phase. The boundary condition at x = a\ directly influences the boundary 
condition at x — b\, in a way mathematically expressed by 

?*(*,&*-) = g^t-^,^ (2.8) 

where Li is the length of the link, k{ is the speed of forward wave propagation. See 
Figure [5] for an illustration. 

ii) the shock reaches the entrance, i.e. x*(t) = a\. In this case, the current link is dominated 
by congestion, and the boundary conditions at x — b{ directly affect the boundary 
condition at x = 

ft (t, a,+) = qAt- — , b-] (2.9) 



where W{ is the speed of backward wave propagation. 




Figure 5: Example of latent separating shocks. Case i): the separating shock reaches the right 
(downstream) boundary. Case ii): the separating shock reaches the left (upstream) boundary. 

In either case, the demand-supply functions will involve boundary flows ( |2.9| ), ( |2.8[ ), which 
are exogenous. The next key step towards the Link-based Kinematic Wave Model is to detect 
these two extreme cases using only boundary flows ^(-, a^), ^(-, hi). This can be done using 
the Lax-Hopf formula [H E3 EH B3 S2] • 



2.4 The Lax-Hopf formula approach for detecting latent shock 

We will present, in this article, another solution method known as the Lax-Hopf formula. 
Initially introduced in [321 [33], then extended in [TJ El 157] . and applied to traffic theory 
in [HI HU [15], the Lax-Hopf formula provides a new characterization of the solution to the 



9 



hyperbolic conservation law and Hamilton-Jacobi equation. The Lax formula is derived from 
the characteristics equations associated with the Hamilton-Jacobi PDE, which arise in the 
classical calculus of variations and in mechanics, see [T7] for a complete discussion. 
In view of the following conservation law on link 

d d 

let us introduce the function A^(-, •) : [0, +oc) x [a, b] such that 

d d 

faNi(t,x) = -Pi{t,x), Q- t N ^ x ) = fi(Pi(t> x )) ( 2 - 10 ) 

One way of interpreting Ni(t, x) is that it is the cumulative number of vehicles that have 
passed location x by time t. The function N{(t, x) is sometimes referred to as Moskowitz 
function or Newell-curves and has been studied extensively, in [HI El [HJ EJ S2] . A well-known 
property of A^(-, •) is that it satisfies the following Hamilton-Jacobi equation 

^Ni(t, x) - fi (-^-Ni(t, xfj = 0, (t, x) E [0, +oc) x [a,, b t ] (2.11) 

Consider the triangular fundamental diagram given by 

Hp) = i hp pe[0 >P J] 

h[P> \-w t (p- P r m ) p£(pt,pT} ' 

where k{ E 5ft+ is the speed of forward propagation of waves, while W{ E is the speed of 
backward propagation of waves. p\ is the unique density such that fi(p*) = Ci, Ci is the 
flow capacity. Based on this particular fundamental diagram, let us introduce the concave 
transformation of the fundamental diagram fi(-): 

f*(u) = sup {fi(p)~up} = d-p*u, ue[-Wi,ki] (2.13) 
pe[o,pl am ] 

Before we articulate the Lax-Hopf formula, let us define the upstream and downstream 
boundary conditions. 

Definition 2.4. (Boundary conditions for H-J equation ( |2.1lD ) Consider a link I i with 
flow capacity Ci. The upstream boundary condition N up ^(-) : [0, +oo) and the down- 

stream boundary condition Nd OW n,i [0, +oc) — >► are non- decreasing, Lipschitz continuous 
functions with Lipschitz constants Ci, such that N up ^(0) = Nd OW n,i(0) — 0. 

Note that the last condition in the previous definition is to ensure zero initial condition on 
the link. The following theorem (Lax-Hopf formula) has been derived separately in various 
contexts in [HJ US [37] 

Theorem 2.5. (Lax-Hopf formula) Consider the Hamilton-Jacobi equation (2.11]) with flux 
function (2.12fy . Let Qi be a subset of the domain [0, +oc) x [a^, bi] and q(-, •) : Qi —> 
be the value condition which prescribes the value of A^(-, •) on Qi. Then the solution to the 
Hamilton-Jacobi equation \2. 1 1] with value condition q(-, •) is given by 

Ni(t,x) = wf{ci(t-T,x-Tu) + T fi(u)} (2.14) 

such that u E [—w^ ki], T > and (t — T, x — T u) E fV 
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Next, we will instantiate the general Lax-Hopf formul (2.14) to the special case with two 



boundary conditions, that is, the case where Qi = [0, +oc) x {a^, 6^}, and •) is given by 
N up ,i( m ) and Nd OW n,i(') • The next result was given in [32]. 

Proposition 2.6. (Instantiation of Lax-Hopf formula for boundary conditions) 

Given upstream and downstream boundary conditions N uPi i(-), Ndown,i(') as ^ n Definition 2.4 



the solution to the H-J equation (2.11) is given by 



Ni(t, x) 



min < N 



t- 



ki 



N do 



t - 



+ P> am (bi- 



(2.15) 



Remark 2.7. The geometric meaning of simplified expression (2.15) is clear. Since the tri- 



angular fundamental diagram implies only two characteristic speeds, the point (£, x) can be 
the intersection of at most two characteristic lines, one with positive speed k, another with 
negative speed —w. Thus the Lax-Hopf formula reduces to the comparison between these two 
cases. Such insight was provided by Newell 

Recall that the separating shock x*(-) : [0, +oc) — >► [a, b] is a continuous curve in the 
temporal-spatial domain. The next proposition is the key for detecting latent shocks and 
switch of regime variables r^. 

Proposition 2.8. Let A^(-, •) : [0, +oc) x [a^, bj\ be the unique viscosity solution to Hamilton- 



Jacobi equation (2.11) with zero initial condition and appropriate upstream and downstream 
boundary condition N uPi i(-), Ndown,i{')- Then the following statements hold: 



x — a* 



ki 



1. For any (£, x) E [0, +oc) x [a^, bj\, 

x < x*(t) iff N upA 
x > x*(t) iff N up ^ 

2. In particular, for all t >0, 



kj 



— Ndown,i 



Wi 



+ f ?- m (b i -x) (2.16) 
+ ^ m (b l -x) (2.17) 



<{t) = at iff N uP:i (t) > N, 



down A 



Wi 



■ jam T 

+ (4 L i 



x*(t) = h iff N r 



Li 

ki 



(2.18) 
(2.19) 



Proof. By standard PDE theory [17], the quantity shown in the right hand side of (2.14) is 
minimized along a path p in the temporal-spatial domain if and only if that path is precisely 
a characteristic line connecting the point (£, x) to the boundary. In addition, fix any t, the 
point x is to the left (right) of the separating shock if and only if (£, x) is connected to the 
left (right) boundary ai (pi) by a characteristic line, see Figure [5j Thus fl2~T6l ), ( |2~T7| follow 
from the above two observations. □ 



3 The DAE system 

The dynamics described via the separating shock and Lax-Hopf formula (Section |2.4| ), together 
with the appropriate Riemann Solver (Section |2.2|), can be explicitly written as a system of 
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differential algebraic equations. The DAE system is a common mathematical form of network 
loading problems arising from analytical dynamic traffic assignment [20j |2T] . The fact that 
the LWR-based network flow models - which are usually viewed and studied as a constrained 
partial differential algebraic equation (PDAE) system [3lJ SS] ~~ can be re-formulated as an 
explicit DAE system, provides new ground for mathematical analysis and computation. 

Let us consider a road network represented by a directed graph (A, V), where A, V denotes 
the set of links and nodes, respectively. It is not restrictive to assume that all the nodes are 
of the type depicted in Figure [2] (notice that one may consider nodes with more complicated 
topology upon acquisition of appropriate Riemann Solvers, which does not change the nature 
of the current DAE formulation). We begin by introducing the following notations for each 
link L E A: 



Qin,i( 
Qout,i ( 

A( 
Si( 

N U p,i ( 
^down,i ( 



the flow of cars entering link 1^ 

the flow of cars exiting link 

the maximum flow allowed to exit the link 

the maximum flow allowed to enter the link 

the cumulative number of cars that have entered link Ii 

the cumulative number of cars that have exited link 



In addition, let the sets of merging nodes and diverging nodes be D, respectively. For each 
v E V, denote by l v and O v the set of incoming links and outgoing links, respectively. For 



the compactness of notation, we will denote the Riemann Solvers discussed in Section 2.2.2 



and |2.2.I] by RS± and RS2, respectively: 

(flout.ii Qinj) — RS\(Di) 5j), 
(flout, ii Qinj) — RS2(Di, Sj), 



We assume that vehicles try to flow into the network from the set of origins V Q C V, at a 
departure rate h Vo (-) for each v Q E V G . For the consistency of notations, and also for the easy 
transition to a Cauchy problem on networks, we assume that each origin v Q is connected to 
an virtual link I Vo E I Vo that is sufficiently long. Vehicles are assumed to be traveling along 
such virtual link before then can enter the network through the origin. However, the actual 
departure rate will be subject to the constraints at vo described by the junction model. Such 
consideration is necessary since when the spillback occurs on the links immediately downstream 
to the origin node, not only the travelers from upstream links but also the travelers departing 
from the origin node will be affected. As a result, the desired departure rate h VQ (t) may differ 
from the actual departure rate. 

Using ( |2.5| )-(2.6) and ( |2.18 )-(2.19), we can easily derive the following DAE system: 



^^down,i(t) Qout,i(t)i 



he A 



Di(t) 



Qout,^ 



ki 



if K, 
if K 



up,i 



up,i 



Li 
hi 



>N, 



down, 



(- 



if N up>j (t) = N d 



,ownj 



t - 



(*) 
(*) 

Li 



if N up j(t) < Ndown,: 



f.er 



, jam T 

+ Pj Lj 



(flout, ii Qinj) — RSl(Di, 5j) 5 
(flout, ii Qinj) = ^(A 5 Sj)i 



h 

he I" 



v e M 
v g V 



(3.20) 
(3.21) 

Ij € O v (3.22) 

(3.23) 
(3.24) 
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4m,j — 



where Li is the length of link 1{. k{, W{ are the forward and backward wave propagation 
speeds, respectively. pj am is the jam density on link I{. 

4 Analysis of the continuous-time model 



The DAE system (3.20)-(3.25) is very convenient for computation. In fact, a direct time- 



discretization of the system coincides with the Link Transmission Model proposed earlier in 



[5T] . On the other hand, due to the numerous discontinuities emerging from (3.21), (3.22) and 



the Riemann Solvers (3.23), (3.24), the DAE formulation is not convenient for mathematical 
analysis. In particular, it is difficult to show solution existence and uniqueness, as well as 
continuous dependence on the initial/boundary conditions, by directly analyzing the DAE 
system. 

This section is dedicated to the analytical properties of the continuous-time solution men- 
tioned above. In order to investigate our model in a framework that allows rigorous mathe- 
matical analysis, we will use a wave- front tracking algorithm as means of approximating the 
solution. 

4.1 A wave- front tracking algorithm on networks and existence of solutions 

A wave-front tracking algorithm consists in a procedure to construct piecewise constant ap- 



proximation of solution to equation ( 1.1 ) on the network. We briefly summarize here the main 
steps, referring to [3j for a detailed description and to [22J for the network case. 

Given a discretization parameter 6 and initial-boundary conditions in BV - the space of 
functions with bounded variation, an approximate solution is constructed in the following way: 

• Approximate initial-boundary conditions by piecewise constant functions and solve the 
Riemann Problems (RPs) at discontinuities of such approximations and at junctions. 
Approximate rarefactions by rarefaction shocks of size 6; 

• Use the piecewise constant solution obtained piecing together the solutions to RPs up to 
the first time of interaction of two travelling waves (shocks) or of a wave with a junction; 

• Then solve a new RP created by the interaction and prolong the solution up to next 
interaction time, and so on. 

To ensure the feasibility of such construction, it is enough to control the number of waves 
and interactions. This is easily done in scalar case on a single link since both the number of 
waves and the BV norm are decreasing in time (see [3J for details). For networks, one has 
to carefully estimate all quantities, due to interaction of waves with junctions. Indeed, when 
a wave interacts with a junction from a link, it may produce new waves in all other links 
connected to the junction. For the network case see [22| [23] . 

Now consider a sequence of approximate solutions p$ constructed by a wave-front tracking 
algorithm. If one can provide estimates on the total variation, then as S tends to zero p$ tends 



to a weak entropic solution of (1.1) on the whole network. Technically one uses compactness 



in BV provided by Helly's Theorem and the weak formulation of equation ( |1.1[ ) (see [3} [22]). 
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Definition 4.1. Fix an approximate wave front tracking solution p§ and a node J. An in- 
coming link Ii is said to have a good datum at J at time t if 

n(t, hi-) 1 

and a bad datum otherwise. Similarly, we say that an outgoing link Ij has a good datum at J 
at t if 

rj(t, aj +) = 

and a bad datum otherwise. 

In [23] a general existence theorem (Theorem 5.1) for the weak solution at a junction was 
provided based on properties of the given Riemann Solver. More precisely, three sufficient 
conditions for existence are given: 

(PI) The solution to Riemann problems depends only on the values of bad data. 



(P2) For waves interacting with a junction, the change in flux variation (i.e. the sum of flux 
jumps over all waves) due to the interaction is bounded by the change in flux through 
the junction and by the flux jump of the interacting wave. 



(P3) Interactions of waves bringing a flux decrease produce a decrease in the flux through the 
junction. 

Proposition 4.2. The Riemann Solvers for merge and diverge junctions in Section \2l^ satisfy 
properties (Pl)-(PS), thus the weak solutions to the Cauchy problem at the junctions exist. 
Furthermore, the global solution to the Cauchy problem on a network exists. 

Proof. The case of merge was already treated in [23], thus we focus on the diverge case. 



From formulae (2.5), (2.6) and Definition 4.1, we have that supply and demands are determined 



only by bad data, thus property (PI) follows. 



From (2.7), we can deduce that the change in the flux through the junction is proportional 
to the flux jump of the interacting wave. Moreover, the flux jumps in all roads, after the 
interaction, is proportional to the change in the flux through the junction. Thus (P2) follows. 
Finally, a wave bringing a decrease of flux produces a decrease of supply-demand functions, 
thus, again by ( |2.7| ), a decrease in the flux through the junction. So (P3) is granted. Thus 
Theorem 5.1 of [23J implies that a weak solution exists at the junction. 

To show global existence of a weak solution, one only needs to notice that the key feature 
of conservation law models is the finite speed of propagation. Thus if the local existence is 
guaranteed, so will be the global existence. □ 



Remark 4.3. The existence result in Proposition ^ • S\ also applies to boundary value problems 
on networks. The boundary conditions can be easily transformed into an initial condition by 
introducing virtual links, as discussed in Section [3| 

4.2 Well-posedness and uniqueness of solutions 



To prove uniqueness and well-posedness for ( 1.1 ) one may rely on techniques such as Kruzhov's 



entropies (see [H]), Bressan-Liu-Yang functional (see [3J) or generalized tangent vectors (see 
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[22]). For networks, it happens that generalized tangent vectors appear to be the most con- 
venient tool. 

Given a piecewise constant function / : [a, 6] —> i?, a tangent vector is given by shifts of 
the discontinuity points of /. More precisely, let us indicate by U the discontinuities of /, 
where to = a < t\ < • • • < tjy < £/v+i — b, and by fi the value attained by / on [U-i^U]. 
Then a tangent vector to / is given by £ = (£1, . . . , £/v) E R N . For each e, the corresponding 
perturbation of / is given by: 

fe(x) = fi x E [U-i + e£i-i,U + £&[ 

for i = 1, . . . , N + 1, where we set £o £/v+i 0- 
The norm of £ is defined as: 

N 



2=1 

in other words the tangent vector norm is the sum of norms of ^ multiplied by the size of the 
shifted jumps. 

To prove uniqueness and well-posedness one proceeds as follows. Given piecewise constant 
initial-boundary conditions, one considers tangent vectors to them. Then one proves that the 
norms of tangent vectors remain bounded along approximate wave- front tracking solutions. 
This implies that L 1 distances are uniformly bounded, see [22] for details. 
The key fact in such procedure is that for interactions of waves with junctions the norm 
of tangent vectors changes in the same way as the jumps in the fluxes. More precisely, let 
us consider a junction J and assume the interacting wave (p^~, pf) is coming from road I{. 
Indicate by Aqi the jump in the flux and by ^ the shift. Also indicate by Aqj the jump in the 
flux on some other road Ij resulting from the interaction and by £j the corresponding shift. 
Then we have the following lemma ([22J) 

Lemma 4.4. If the wave on I{ interacts with J without producing waves in the same road 
then the shift £j produced on Ij, as a result of the shift ^ on I{, satisfies 

m-Pj) - ^(Pt-PT) (4-26) 
where pj, pj are the states at J on road Ij before and after the interaction, respectively. 
Proof. By Rankine-Hugoniot condition, the speed of wave (pf, p^) is 

f(pf)-f(p-) 
pi - p7 

If the wave (p^ , pf) is shifted by £j, then the interaction of this wave with J is shifted in time 
by 

Pi -Pi _ c Pi -Pi 



W)-ttp7) A * 

Since the shift in time taking place on Ij must be the same, we have 

r Pi-P7 „ p pI-pJ 



□ 
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Notice that the quantity (p+ — p^~) | is precisely the L 1 -distance of two initial conditions 
on Ii with and without the shift ^ at the discontinuity (p^, p+), respectively. Same situation 
applies to other roads Ij. Rewrite ( 4.26[ ) as 



Sj(pj-Pj)\ = ^\^Pt-PT)\ (4-27) 

Therefore, to control the norm of tangent vectors one has to check that the multiplication by 
factors i^|4 remains uniformly bounded, regardless of the number of interactions which may 
occur with the junction. (This bound may indeed fail for general junctions, see [22J, Section 
5.4, Proposition 5.4.4). 

We are ready to prove the following: 



Proposition 4.5. For Merge and Diverge junctions in Section 2.2, the increase in tangent 



vector norms as in (4-27) is uniformly bounded. 



Proof. The case of Merge junction was already treated in [23] . so we focus on Diverge junction. 
Consider a wave interacting from road i and the one produced on road j, where i, j = 1, 2, 3. 



By (2.7), the variations of the fluxes A^, i = 1, 2, 3 satisfy 

A<?2 = ai j2 Agi, Aq 3 = a^ 3 Aq 1 
consequently, the multiplication factors satisfy 

Ai A 2 A 3 

A^ ' A^ ai ' 2 ' AT = ai ' s; 

Ai _ _J_ A 2 _ A3 _ 0^3 

A 2 0*1,2' A 2 ' A 2 ai ;2 ' 

Ai 1 A 2 a^2 A 3 



A 3 ai j3 ' A 3 ai >3 ' A 3 



1. 



Then the multiplication factors 7^ of tangent vectors are given by the following table (rows 
i and columns j): 

1 ai,2 «i,3 

CK1,2 CK1,2 (4.28) 

1 Ql,2 2 

«1,3 CKl,3 

Since 

7ij7jfc = Tiife? h h k = 1, 2, 3 

we have that the norm of tangent vecotors are uniformly bounded independently from the 
number of interactions with the junction. □ 
With the uniform bound on the norm of tangent vectors, one proceeds in a standard way 
as in [31 [22] to show the well-posedness of Cauchy problem with any BV initial/boundary 
data. The key is to construct a one-parameter curve joining two initial data, then estimate 
the length of the curve by integrating the length of the tangent vector. 

Remark 4.6. The argument put forward in this section can be extended to treat networks 
with origin- destinations. The idea is to replace distribution matrix A with piecewise constant 
matrix A(-), and the established results will apply. However, in the process of passing to the 
limit, the elements of A could be arbitrarily small, causing the norm of tangent vectors to blow 
up. A detailed analysis of this situation is left as future research. 
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5 Numerical example 



In the numerical example, we consider a simple network shown in Figure [6] consisting of seven 
arcs I\ — I7 and four nodes a, 6, c, d. The parameters associated with each link are displayed 
in Table 1 




Figure 6: A network example consisting of seven arcs and four nodes 
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Table 1: Arc parameters for the small network. p^ am : jam density; h{\ forward wave speed; wf. 
backward wave speed; Cf. flow capacity; length. 



We consider a time horizon [0, 5] spanning five hours, with a time step 5t = 0.05 hour. We 
assume vehicles enter the network from link I\ during a time interval [0.5, 3.75]. The inflow 
on each time interval is a random variable that is uniformly distributed between and the 
flow capacity C\. A plot of the inflow is shown in Figure [7[ 

The network only contains two types of junctions: the ones with one incoming link and two 
outgoing links, and the ones with two incoming links and one outgoing link. For simplicity, 
we set both distribution parameters at diverge junctions and the right-of-way parameters at 
merge junctions to be 1/2. 

As indicated by Table 1, arc I4 is a bottleneck. In addition, the merge node c also suggests 
the existence of a potential bottleneck since the capacity of downstream link Iq is less than 
the sum of capacities of I3 and I4. In order to have a clear visualization of the backward 
propagation of waves and the separating shocks on each link, we use the boundary data 
Qout,i, Qin,ii i = 1, . . . , 7 obtained from the solution to construct the Moskowtiz function via 
Hamilton- Jacobi equation (2.11) and Lax-Hopf formula. For the H-J equation, the separating 
shock no longer represents discontinuity, rather, it is displayed as a 'kink' (discontinuity in 
the first derivative). The Moskowitz function for links /i,/2 5 ^3 and I4 are shown in Figure pi 
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Figure 7: Inflow profile on link I\. 



[9j[T0]and[TT 



From these Moskowitz functions, we can clearly see two types of regimes: the free-flow 
phase and the congested phase. The mutual boundary of these two regimes is the separating 
shock. One can observe that the congestion initiates from node c (downstream boundaries of 
link ^3 and I4) and propagates backward along links 1% and I4, then further to the upstream 
links I\ and I2. 

We also observe that on links ^3, ^4, there exist latent shocks on the upstream boundaries 
as and for some period, indicating the existence of spillback. To have the most direct 
visualization of the spillback phenomena, we put the Moskowitz functions of /1, ^3 (respectively 



I2, I a) side by side, and viewed from above. These are displayed in Figure 12 (13), which clearly 
shows the formation, propagation and dissipation of congestion on a road network. 
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Figure 8: Moskowitz function for link I\. Figure 9: Moskowitz function for link I2. 
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Figure 10: Moskowitz function for link 1%. Figure 11: Moskowitz function for link I4. 



Link I 



LinkT 




Figure 12: Moskowitz functions on links I\ and Is- The congestion on link ^3 spills onto link 
h. 

6 Conclusion and future research 



This paper presents a continuous-time traffic flow model on vehicular network. We use flow 
q and regime r as state variables, and implicitly solve the dynamics inside each link using 
the notion of separating shock and the method known as Lax-Hopf formula. On top of this 
formulation, a DAE system can be formulated explicitly describing the temporal evolution of 
the system states. The DAE system is the continuous-time version of the Link Transmission 
Model [51]. Analytical results of the system such as solution existence, uniqueness and well- 
posedness of Cauchy problem are established. 

As our next step towards an analytical framework of network flow models and DTA models, 
we will consider networks with Origin-Destination and the users have prescribed route choices. 
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Figure 13: Moskowitz functions on links I2 and I4. The congestion on link I4 spills onto link 

h. 



As mentioned in Remark 4.6, such model poses new challenge in the estimation of norms of 



tangent vectors in the front tracking approximation. The DAE system (3.20)-(3.25) also needs 
to be extended to the O-D network. 
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